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Abstract. New models of rotating and non-rotating stars are computed for initial masses between 25 and 120 
Mq and for metallicities Z — 0.004, 0.008, 0.020 and 0.040 with the aim of reexamining the wind contribution 
of Wolf-Rayet (WR) stars to the ^^F enrichment of the interstellar medium. Models with an initial rotation 
velocity Vi — 300 km s~^ are found to globally eject less '^^F than the non-rotating models. We compare our new 
predictions with those of Meynet & Arnould (2000), and demonstrate that the ^®F yields are very sensitive to 
the still uncertain ^^F (a , p) '^'^Ne rate and to the adopted mass loss rates. Using the recommended mass loss rate 
values that take into account the clumping of the WR wind and the NACRE reaction rates when available, we 
obtain WR ^^F yields that are significantly lower than predicted by Meynet & Arnould (2000), and that would 
make WR stars non-important contributors to the galactic ^^F budget. In view, however, of the large nuclear and 
mass loss rate uncertainties, we consider that the question of the WR contribution to the galactic ^^F remains 
quite largely open. 
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1. Introduction 

The solar system has for long been the only location in 
the Universe with a known fluorine (^^F) content. The 
very origin of this ^^F has been a major long-standing 
nucleosynthetic puzzle, in spite of the fact that it has the 
lowest solar abundance among the nuclides ranging from 
carbon to calcium. 

Since the beginning of the nineties, the situation has 
changed drastically, both theoretically and observation- 
ally. The first quantitative prediction that "'^^F could be 
produced thermonuclearly at a level compatible with the 
solar amount has been made by Goriely et al. (1989). They 
identify a mode of production of ^^F that could develop 
in He-burning conditions following the reaction chains 

{l3+)''0{p,a)^'N{a,j)^^F 

/ \ 
i4iV(a,7)i«F ^ (n,p)i80(p,a)i5iV(a,7)i9F(a,]5)22iVe. 



(n,a)i57V(a,7)i^F 



In this scenario, the synthesis of ^^F requires the availabil- 
ity of neutrons and/or protons. They are mainly produced 
by the reactions (a , n) ^^^O and (n , p) . 
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The viability of the proposed He-burning ^^F produc- 
tion scenario has been demonstrated in the framework of 
detailed models for aymptotic giant branch (AGB) stars 
(Forestini et al. 1992, Mowlavi et al. 1998, Lugaro et al. 
2004), as well as of massive stars evolving through the 
Wolf-Rayet (WR) stage (Meynet & Arnould 1993, 2000 
hereafter Paper I; Stancliffe et al. .2005i) . Massive stars that 
do not experience the WR phase are expected to produce 
instead an insignificant amount of ^^F during their hy- 
drostatic evolution (Meynet & Arnould 1993, Woosley & 
Weaver 1995, Limongi & Chieffi 2003). Note that the Type 
H supernova explosions of massive stars have also been 
claimed to be responsible for a ^^F production through /z- 
and T-neutrino spallation on ^"Ne (e.g. Woosley & Weaver 
1995). This production is highly uncertain, as it is very 
sensitive to the poorly known neutrino energy spectra. It 
will not be discussed here. 

The most recent solar and meteoritic ^®F values are 
provided by Fodders (2003) and Asplund et al. (2005). 
Observational efforts to determine ^^F abundances out- 
side the solar system have been largely triggered by the 
early theoretical predictions of Goriely et al. (1989). In 
a companion paper to the one of Forestini et al. (1992), 
Jorissen et al. (1992) provide the first ^^F abundances 
measured in stars other than the Sun. They analyse a 
number of s-processed enriched galactic MS, S, and N- 
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type giants having a near-solar metallicity. These observa- 
tions demonstrate that AGB stars are fluorine producers, 
nicely confirming the initial predictions by Goriely et al. 
(1989) and Forestini et al. (1992). The derived ^^F abun- 
dances correlate with the carbon and s-nuclide ones, a 
pattern that AGB models can also account for, as first 
shown by Mowlavi et al. (1998). Very recently, ^^F has 
also been detected in a sample of H-rich as well as H- 
deficient PG1159-type hot post- AGB stars (Werner et al. 
I2()()4|l . The H-rich stars, which are not especially C-rich, 
show solar-like ^®F abundances. This is in line with the 
conclusion by Jorissen et al. (1992) concerning the C-F 
correlation in AGB stars. In contrast, the H-poor stars 
exhibit very large F overabundances ranging from 10 to 
250 times solar. Werner et al. (2004) suggest that this 
abundance pattern might be explained by the operation 
of a late post- AGB shell flash (see Herwig et al. 1999). As 
a complement, a F/H abundance ratio of 4.5 x 10^^ has 
also been determined in a planetary nebula (Liu 1998). 
This value again agrees with the Jorissen et al. (1992) 
correlation between F and C. 

Fluorine data also exist for stars whose surface, in con- 
trast to the (post-)AGB stars, is not expected to be con- 
taminated with the products of their in-situ nucleosyn- 
thesis. This concerns a few near-solar metallicity K and 
M giants also analysed by Jorissen et al. (1992). In ad- 
dition, Cunha et al. (2003) have studied a sample of red 
giant stars in the Large Magellanic Cloud, as well as in the 
atypical galactic globular cluster lo Centauri. These low- 
metallicity giants exhibit sub-solar fluorine abundances 
(A(F)i e [3; 4.19]) , while both the K - M field giants and 
the three K - M pre-main sequence low-mass stars of the 
Orion nebula cluster exhibit nearly solar fluorine (A(F) 
~ 4.55) (Cunha & Smith 2005). This effect of metaUic- 
ity is in agreement with the general behaviour of [F/0] 
exhibited by stars in different evolutionary phases within 
the galactic disk. Such observations of non-contaminated 
stars spanning a range of metallicities are mandatory if 
one wants to build a model for the evolution of the ^^F con- 
tent of the Galaxy and of other stellar systems (Renda 
et al. 2004). Finally, fluorine has also been observed in 
various interstellar locations (Federman et al. 2004, and 
references therein) in an attempt to constrain the ^^F nu- 
cleosynthesis models. 

The aim of this paper is to revisit the predictions of 
Paper I of the ^^F yields of WR stars on grounds of new 
models for different masses and metallicities, and to pro- 
vide the first predictions of the ^^F production by rotating 
WR stars. 

2. Input physics 

The models used here are computed with the Geneva 
stellar evolution code from the Zero Age Main Sequence 
(ZAMS) up to the end of the He-burning (HeB) phase and 
are listed in Tab. ^ The physical ingredients, structural 

1 A(X) = log[n(X)/n(H)]-|-12 



predictions, and comparisons to observations are discussed 
at length by Meynet & Maeder (2003, Figs. 9-12 and Tabs. 
1-2 ;|2nnil Figs. 7-9 and Tables 1 and 3). The main points 
of relevance to the ^^F synthesis are as follows: 

(1) The initial compositions for the different metallici- 
ties are selected as in Palacios et al. (2005). The initial 
^^F mass fraction Xiq = Xiq{Z,{)) at metallicity Z is de- 
rived from the simple scaling Arig(Z, 0) — (Z/0.02)Xi9o, 
where 0.02 is the adopted metallicity for the Sun, and 
XiQQ = 4.1 X 10~^ is the solar ^^F mass fraction accord- 
ing to Grevesse & Noels l)1993|l : 

(2) The effect of rotation on the mass loss rate M is taken 
into account as in Maeder & Meynet 1)2001(1 . As reference 
M, we adopt for the pre-WR stages the values proposed 
by Vink et al. (200131201111) who account for the occurrence 
of bi-stability limits which affect the wind properties and 
mass loss rates. Outside the domain covered by these au- 
thors, the rates from de Jager et al. H1988|l are selected. 
As the empirical M values are derived from stars with 
a variety of rotation velocities, and as M decreases with 
these velocities, a reduction factor to the empirical rates of 

0. 85 (Maeder & Meynet 120011 is introduced for the non- 
rotating models. During the WR phase, we use the M pre- 
scriptions of Nugis & Lamers H2000|l . These rates, which 
account for the clumping of the winds, are 2 to 3 times 
smaller than the ones used in previous non-rotating 'en- 
hanced mass loss rate' stellar models presented in Paper 

1. Note that wind anisotropics induced by rotation are ne- 
glected. These anisotropics are indeed shown to be very 
small for the initial velocity Vi — 300 km s^^ (Meynet 
& Maeder 12003 selected in this work (see point 6 below). 
This would not be true for higher initial velocities (Maeder 
!»•)): 

(3) During the pre-WR phases, it is assumed that the 
mass loss rates have a metallicity dependence given by 
M{Z) = (£/0_02)1/2m(0.02) (Kudritzki & Puis IWiHl 
Vink et al. 1200 l|l . In contrast, no metallicity dependence 
is introduced during the WR stage; 

(4) All the models are computed with moderate core 
overshoot. The distance of overshoot is taken equal to 
d = aHp, where Hp is the pressure scale height at the 
Schwarzschild boundary and a — 0.1. This value of a is 
twice as small as the value used in the models of Paper I; 

(5) The transport of the nuclides and of the angular mo- 
mentum is described as in Maeder & Meynet pM|l and 
Meynet & Maeder p]?)^ ; 

(6) All the considered stars are assumed to rotate on the 
ZAMS at an initial rate vi = 300 km s'^ . For Z = 0.02, 
this value leads to time averaged equatorial velocities on 
the Main Sequence (MS) well in the observed range (be- 
tween 200 and 250 km s^^ ); 

(7) The reaction rates adopted in Paper I have been 
updated by the use of the NACRE data (Angulo et 
al. 11999(1 . when available. The rates of the reactions 
entering the chain displayed in Sect. ^ ^^^id that are 
not considered in NACRE arc taken from the fol- 
lowing references: ^-^N (n , p) ^-^C (Brehm et al. HM^ . 
i^F(n,p)i8o (REACLIB, Thielemann et al. HW!^ . 
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Table 1. Values of p™g"'^(Mi, Z) in units of lO^^Mo for the rotating and non-rotating models also displayed in Fig.O 



M(Mq) 


Wi (km s ^ ) 


Z = 0.04 


prr (10-" Me) 

Z = 0.02 Z = 0.008 


Z = 0.004 


120 


300 


5.38 


-24.13 18.6 


29.1 


120 





51.86 


110.05 




85 


300 


28.41 






60 


300 


27.97 


19.4 -5.18 


-1.81 


60 





103.5 


22.96 




40 


300 


27.67 


-2.32 


-0.757 


30 


300 




-1.49 


-0.284 


25 


300 


15.17 


-6.03 




25 





-1.6 


-0.201 





i^F (n , a) (Caughlan & FowlerlMi hereafter CF88), 
and i9F(a,p)22Ne (CF88). 

3. The WR production of ^^F 

Figure ^ shows the abundances of various nuclides inside 
the non- rotating 60 M0, Z = 0.02 model at three stages 
during the core He-burning (HeB) phase associated with 
(a) the maximum mass fraction Xig^c of fluorine at the 
centre, (b) the first time the mass fraction Xig.s of fluo- 
rine at the surface exceeds 10~^, and (c) the maLximum 
mass fraction X19 s of fluorine at the surface. Panel 
(d) refers to the beginning of the C-burning phase, 
when there is no fluorine left in the core. We focus 
on the HeB stage since ^^F is produced in the convective 
core at the beginning of that phase mainly through the 
first chain of transformation of ^^N into ^^F displayed in 
Sect, n In only about 43 300 yrs, while the central He 
mass fraction Xi^c drops from 0.98 to about 0.90, the 
^^F mass fraction at the centre Xig.c is seen to increase 
from 4.4 x 10~^° to its maximum value 1.6 x 10~^ (see 
panel (a)). At this time, the surface ^^F mass fraction 
Xig^s is still very small (3.6 x 10"^"). As evolution pro- 
ceeds, the convective core retreats in mass, leaving in its 
wake ^^F-enriched layers that will eventually appear at 
the surface when the star enters the WC phase (see pan- 
els (b) and (c) ) . The strong wind during this phase finally 
allows the ejection of these fluorine-rich layers in the in- 
terstellar medium (hereafter ISM). Panel (b) depicts the 
situation when Xig,s just exceeds 10~^. The time elapsed 
between panels (a) and (b) is of the order of 148 000 
yrs, corresponding to a X4^c decrease from 0.90 to 0.47. 
During this period, ^^F in the convective core is partly de- 
stroyed by i9F(a,p)22Ne. In panel (b), Xig^c = 3.7x 10-^ 
which is more than 4 times lower than its maximum value. 
Concomitantly, mass loss starts exposing the most ^^F-rich 
layers at the surface (panel (c)). From this stage on, the 
synthesised "'^^F starts enriching the ISM until the outer 
layers are severely -^^F-depleted, as it is the case at the 
beginning of core C-burning (panel (d)). 

A generic sequence summarising the above can be 
used to describe the evolution of the central and sur- 



face ^^F mass fractions in all the models presented here. 
Fluorine is produced in the convective core and its mass 
fraction rapidly increases at the beginning of the He- 
burning phase. The regions above the convective core 
are not ^^F -rich at this stage. As the evolution pro- 
ceeds, the central temperature becomes high enough for 
^^F (a , p) ^^Ne to be efficient, leading to a decrease of the 
^^F mass fraction in the convective core. On the other 
hand, the convective core retreats in mass as a result of 
strong mass loss. This allows part of the fluorine left be- 
hind by the retreating core to escape destruction. This flu- 
orine can then be exposed at the surface by the stripping 
of the outer layers by stellar winds. Thus, at the surface, 
the ^^F mass fraction increases, reaches a maximum, and 
then decreases as deeper layers are revealed. 

From the above, it appears that WR stars can be 
^^F contributors to the ISM if at least : 

(1) The star can enter the WC phase at a sufficiently 
early phase of core He-burning, so that the ^^F -enriched 
shells of the He-core, which coincide with those rich in 
carbon and oxygen, can appear at the stellar surface be- 
fore fluorine is too much depleted (see Sect. 3). 

(2) The mass loss at the beginning of the WC phase is 
high enough for removing efficiently the ^^F-rich layers. If 
mass loss is weak at this stage, part of the ^^F has time to 
be converted into ^^Ne when deeper layers are uncovered. 
The fulfilment of these two conditions depends on initial 
mass, metallicity and rotation velocity, as discussed below. 

The ^^F ISM enrichment efficiency of a star with initial 
mass Mi and metallicity Z, undergoing a wind phase of 
duration T(A/i, Z) is conveniently evaluated in terms of its 
net i^F yield p^^'"^ defined as 

prt''iM„Z)= (1) 

/ [Xig,s{Mi:Z,t) - Xig,s{Z,0)]M(Mi,Z,t)dt. 

JO 

This quantity is a classical input in galactic chemical evo- 
lution models. However, it does not represent the total 
yields of ^^F , since more fluorine could be produced dur- 
ing the supernova explosion. Note also that negative p™'^ 



4 



A. Palacios, M. Arnould & G. Meynet: WR production of ^^F 





-2 
-4 



GO 

° -6 
-8 
-10 



-| 1 1 1 1 1 1 — 7TTT 



-160 



•H 

14N 



15 



N 



(a) 



191 



I I I I I I 



T I , 



I I I 



20 25 
Mr (Me) 





-2 
-4 
-6 
-8 
-10 



— 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 


— 1 — 1 — r— 1 — 

' 4He J 




- 12C : 

- 160 




'^N 1 

19p - 






~ .... 1 .... 1 ... . 


1 . . . . 



30 17 



18 19 20 
M, (Me) 



21 



X 



-4 - 



_ 1 1 1 1 1 


, , 1 , , 


' 4HeJ 












--^ 160 J 










igp 




. . : ! 1 


: : ! , , 


, 1 , , - 


10 


15 


20 




K (Mo) 





d -8 
-10 




° -6 
-8 
-10 



Fig. 1. Mass fractions of the CNO species and of ^^F versus mass Mr (in Mq) inside a 60 Mq with metaUicity Z = 0.02 
and Vi = km s~^ at three stages during the core He-burning phase (panels a~c), and at the beginning of core C- 
burning (panel d). Panel (a): Xiq reaches its mEiximum value at the centre; (b): Xiq starts to exceed about 
10~^ at the surface; (c): Xiq reaches its mELximum value at the surface; (d): There is no more fluorine 
in the core at the beginning of the C-burning phase. 



values are obtained when the ejected material contains less 
-•^^F than originally present in the star. The yields for the 
wind phase derived for both our rotating and non-rotating 
models are displayed in Table in units of 10~^ M©. 

4. Sensitivity of ^^F production to initial 
conditions 

4.1. Effect of rotation 

Our calculations indicate that rotation has a very lim- 
ited impact on the ^^F production in the stellar cores. 
This largely results from the fact that rotation hardly af- 
fects the temperatures and densities in the core, as well as 
from the inability of rotation to transport ^"^N from the H- 
burning shell into the He core, which would significantly 
boost the ^^F production. The insensitivity of the ^^F con- 
tent to rotation for the 60 M©, Z = 0.02 model is seen 
from a comparison of Figs. and El Let us note that we 
cannot exclude at this point the possibility of a ^'*N mix- 
ing into the core for rotation velocities and metallicities 
different from the ones considered here. 

The limited role played by rotation on the ^^F content 
of the core does not preclude changes in the ^^F yields. 



as shown in Tableland Fig. O In fact, p^™'^ is found in 
all the considered cases to decrease when rotation is in- 
cluded, except in the 25 M© case at Z = 0.04. The mag- 
nitude of this reduction shows a high sensitivity to mass 
and metaUicity. The pj™'^ lowering with rotation can be 
explained as follows. Rotation and associated transport 
processes favour an early entrance of the stars into their 
WR phase (Fliegner & Langer 19951 Meynet & Maeder 
ElOS). Consequently, the period of high mass loss rate is 
lengthened, with the result that the rotating models enter 
their WC phase with, in general, a lower mass. At this 
stage, the mass loss rate scales with the actual luminosity 
of the star (Nugis & Lamers l^OOO|l . and thus with its actual 
mass, since the WC stars obey a mass-luminosity relation 
(Schaerer & Maeder [l992|l . The mass of the ^^F-enriched 
material ejected into the ISM during the WC phase is 
thus lowered as a result of rotation. This p™'^ reduction 
appears to be especially limited for the 60 M© model at 
Z = 0.02. This is due to a balance between different ef- 
fects acting in this particular case (see also Meynet & 
Maeder 120031 . The non-rotating model enters the WR 
phase after a short Luminous Blue Variable (LBV) phase 
characterised by a very high M, while the WR phase with 
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Fig. 2. Same as Fig.[TI but for the Z = 0.02 60 Mq model with vi = 300 km s^^. 



rotation already starts during the MS, so that the LBV 
phase is skipped. The short LBV phase experienced by 
the non-rotating star compensates for its later entry into 
the WR phase. As a net result, the rotating and non- 
rotating 60 Mq models enter the WC phase with about 
the same mass. As far as the 25 Mq , Z = 0.04 case is con- 
cerned, the rotating star ejects some amount of ^^F while 
its non-rotating counterpart does not. This derives from 
the rotating model entering the WC phase, which is not 
the case in absence of rotation. 

In conclusion, mixing induced by rotation does not 
affect the amount of ^^F synthesised in the central re- 
gions of the star at the beginning of the core He-burning 
phase, at least for the range of metallicities explored here. 
This might not be true any more at very low metallicity, 
in which case shear mixing appears to be more efficient 
(Meynet & Maeder EUU^ . Through structural effects (ex- 
tent of convective cores, mass- loss enhancement), rotation 
may however affect the overall evolution of stars during 
the WR phase in such a way as to modify the pY™'^ val- 
ues. 

4.2. EfFect of mass and metallicity 

Stars with different metallicities are likely to produce dif- 
ferent amounts of ^^F. This is confirmed in Fig.Qlfor a ro- 
tating 60 M0 star. As in the non-rotating models, ^19. c 



increases with metallicity (see Paper I). This relates di- 
rectly to the enhanced production of ^"^N by the CNO 
cycles during the central H-burning phase. The larger 
amount of ^'^N available at the beginning of the HeB phase 
allows more ^^F to be synthesised. The very limited change 
of the central densities and temperatures with metallic- 
ity does not affect this conclusion. Increased metallicities 
also favour higher Xig^s values as a result of larger mass 
loss rates {M oc Z^^^) which allow ^^F-enriched layers to 
be exposed at the stellar surface before the eventual par- 
tial transformation of ^^F by ^^F (a , p) ^^Ne . In contrast, 
some ^^F destruction may be unavoidable at lower metal- 
licities. 

These general considerations help interpreting the 
trends of the -'^^F yields with mass and metallicity shown 
in Fig. 13 Let us just focus here on the models with ro- 
tation. For Z < 0.02, < for M-, < 60 Mq. The 
minimum mass for a rotating star to enter the WC phase 
is Afwc ~ 25 and 52 Mq at Z = 0.008 and 0.004. Below 
these limits, the stars cannot contribute to the ^^F enrich- 
ment of the ISM through winds. In the Mwc < M < 60 
Mq range, ^i^^ already dropped below its initial 
value by the time the He-burning products appear at the 
surface, leading to negative p™'^ values. In contrast, the 
Mi > 60 Mq stars contribute to the ^^F enrichment of 
the ISM because enough mass is lost in these cases for 
uncovering the He core before the destruction of ^^F. For 
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the 120 Mq model, vX^'^ is larger at Z = 0.004 than at 
0.008 because the lower- Z star enters the WC phase with 
a higher mass (48 instead of 31 M© for Z increasing from 
0.004 to 0.008), thus losing more mass at this stage. 

The trend discussed above is reverted at Z — 0.02. 
In this case, p^g"'^ grows from negative values for the 25 
M0 star (at Z = 0.02, this star only marginally enters the 
WR phase, and does not contribute to the ISM ^^F en- 
richment) to a positive value for the 60 Mq model. It then 
decreases again to a largely negative value at 120 Mq. For 
this particular case, the model enters the WR phase al- 
ready during the MS and loses more than 80 % of its mass 
prior to the HeB phase. These ejecta are not ^^F -enriched. 
When the wind becomes ^^F-rich at the beginning of the 
WC phase, the star has only a small remaining mass and 
M is consequently small from this point on, allowing only 
a limited amount of ^^F-rich material to be ejected. This 
results in a negative pYg^'^ value. 

At Z = 0.04, pYq^'^ is positive for all models, and is 
seen to be rather insensitive to stellar mass. This results 
from a subtle balance between the dependence with mass 
of XiQ^s and the amount of ejected ^^F-enrichcd material. 

4.3. Uncertainties in the ^^F yields from ^^F(a,p)^^Ne 

The i^F {a , p) ^'^Ne reaction is the main F destruction 
channel in the considered stars (see Sect.^. The large un- 
certainties on its rate of course concur with the problems 
of purely astrophysical nature to affect the reliability of 
the predicted contribution of WR stars to the ^^F bud- 
get of the Galaxy. The status of our present knowledge 
of the ^^F {a , p) ^^Ne rate has been discussed recently by 
Lugaro et al. (2004) and StancHffe et al. (2005). The rate 
is still poorly known, the uncertainties increasing dramat- 
ically with decreasing temperatures. The rate they recom- 
mend is more than one order of magnitude smaller than 
the CF88 rate used in Paper I and in the present work. 
In view of this, p^g"^ ^as been calculated for the 60 Mq, 
Z = 0.02 rotating model decreasing the CF88 rate by a 
factor of 10. As a result, the yield is found to increase by 
more than a factor of 3, as shown by the asterisk in Fig. 13 
The yield could still be increased slightly by decreasing 
the rate further. There appears to be room for this, as the 
lower limit of the rate proposed by Stancliffe et al. (2005) 
is more than 14 orders of magnitude smaller than their 
recommended rate at T = 2 x 10* K! They however come 
to the conclusion that the corresponding yields are only 
increased by at most 10 % if this extremely small lower 
limit is adopted instead of their recommended rate. This 
is not surprising, as the yields are essentially 'frozen' as 
soon as the ^^F lifetime against a-captures becomes longer 
than its residence time in He-burning zones. This situa- 
tion is encountered if the CF88 rate for ^^F {a , p) ^^Ne 
is divided by 2000. Thus, for rates below this limit, p™'^ 
becomes independent of the ^^F a-capture rates, and is 
just some percents higher than those displayed in Tabled 
and Fig. |3 If ^^F (a , p) ^^Ne becomes small enough, one 
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Fig. 5. Same as Fig.Ofor the present non-rotating models 
(open symbols) and from those of Paper I (filled symbols). 
Data obtained for Z = 0.02 and 0.04 are represented by 
squares and triangles. 



might wonder about other '^^F destruction channels, and 
in particular about the precise role of its radiative neu- 
tron captures. In the 60 Mq rotating star at Z = 0.02, 
the mass fraction of neutrons is non-negligible only at the 
very centre of the star, but decreases rapidly by several 
orders of magnitude further out in the convective core. 
Neutron captures are thus not expected to be responsible 
for a significant destruction of ^^F, even if the a-capture 
channel has a reduced efficiency. 

5. Comparison with Paper I 

Figure is similar to Fig. O and presents the yields at 
Z = 0.02 and 0.04 obtained for the non-rotating mod- 
els of Paper I and of the present work. The main differ- 
ences between the two sets of computations lie in the nu- 
clear reaction rates, the amplitude of the core overshoot, 
and the mass loss prescriptions. For the relevant nuclear 
reactions, the rates adopted here are equal to, or differ 
only marginally from the ones used in Paper I, except for 
i5N(a,7)i9F. The rate from CF88 used in Paper I is re- 
placed by the NACRE (Angulo et al. 1999) adopted one, 
which is about 34 times smaller than the CF88 one at 
T — 2 X 10* K. However, this change does not affect the 
^®F production in the non- rotating models. At the begin- 
ning of He burning, the a-captures on ^^N are very effi- 
cient even when adopting the NACRE rate, and ^^N is 
completely transformed into ^^F well before fluorine starts 
turning into ^^Ne. The net amount of ^^F ultimately pro- 
duced and ejected by WR stars is thus controlled by the 
destruction channel ^^F (a , p) ^^Ne rather than by the 
production one. 



Concerning the core overshoot its adopted value is 
lower in the present models than in those of |Paper"l| All 
other things being kept the same, we should expect the fol- 
lowing consequence of this change. Lower core overshoot 
leads to smaller He-cores. Smaller He cores imply that 
more material has to be removed for the core to be un- 
covered, which clearly disfavours "'^^F ejection by stellar 
winds. On the other hand, smaller central temperatures 
result. This does not affect significantly the ^^F produc- 
tion from ^^N, as ^^N has anyway ample time to be totally 
transformed into ^^F. In contrast, the ^^F destruction is 
slowed down, and this tends to make its yields greater. 
Indeed more time is provided to stellar winds to remove 
the outer layers before fluorine is destroyed in the core. 
Reduction of the amount of core overshoot should thus 
have two contrary effects whose relative importance re- 
mains unclear. 

As already mentioned before, the standard models in 
[Paper l| not only differ from the present ones by the core 
overshoot, but by the adopted mass loss prescription as 
well, which also modifles the evolution of the core mass 
and temperature. 

Even if this additional difference prevents us from iso- 
lating clearly the effect of core overshoot, it appears how- 
ever, that the differences with respect to |Paper"T| are dom- 
inated by changes in the mass loss prescription. The pre- 
scription of Sect. El reduces by a factor of about 2 to 3 the 
rates adopted inl Paper"I| It is emphasised above that the 
revised prescription has to be preferred, as it takes into ac- 
count the clumping of the WR wind. In addition, it leads 
to a very good agreement between the observed popula- 
tions of 0-type and WR stars and the predictions relying 
on rotating models (Meynet & Maeder EOni 2003) • 

The impact of a change in the mass loss is intricate, 
and we just try here to identify general trends. As already 
mentioned above, for a given star to contribute to the 
^^F enrichment of the ISM, the mass loss rates have to 
be large enough at the beginning of HeB to uncover the 
core before the ^^F a-particle captures become efficient 
enough. From this, one might expect that pY™'^ increases 
with increasing mass loss rates at the WC phase. However, 
as M scales with the actual mass during this phase, the 
removal efficiency of the "'^^F-rich layers is larger for more 
massive WC stars, which thus need to have lost a relatively 
small mass during the previous evolutionary phases. In 
summary, the WR p^g"'^ yields depend drastically on the 
mass loss prescriptions. They are large only if the mass 
loss rates are high enough for removing most of the outer 
layers at the very beginning of the HeB phase, but low 
enough for the star to keep a relatively high mass when it 
enters the WC phase. 

The general considerations developed above are con- 
firmed by a closer analysis of the non-rotating 60 and 
120 M0 models. As shown in Fig. [S] the former star is 
an extreme illustration of the cases for which the present 
models predict p^g""^ values lower than those reported in 
Paper I, with a reduction larger than a factor of about 
30. Even if, as made plausible by the discussion above. 
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the differences in the predicted yields likely result from 
the combined (and difficult to disentangle) effects of the 
various changes in the ingredients of the two sets of stel- 
lar models, the revised mass loss rates are most probably 
responsible for the new situation encountered for the 60 
Mq model. The Paper I 60 Mq model aX Z = 0.02 indeed 
enters the WC phase with a mass of about 24 Mq and 
Xi,G = 0.79. Only less than 5 Mq remain at the end of 
the evolution. The newly computed model enters the WC 
phase with a mass of about 21 Mq. This mass is close 
to the one computed in Paper I, but is reached later in 
the evolution, at a point where X4.C = 0.47. This delay 
is of course the direct result of the lower mass loss rates 
adopted here in the previous phase. In addition, the final 
mass of the star is 12.4 Mq. In other words, only 8.6 Mq is 
lost during the WC phase, which is less than half of what 
is computed in Paper I. All these effects tend to reduce 

„wind 

Pw ■ 

In contrast to the situation characterising the M < 
60Mq models, the new ^^F yields for the non-rotating 
120Mq model are larger than those of Paper I (Fig.E)). In 
Paper I, the star is predicted to lose so much mass that it 
enters the WC phase with only about 6 Mq . Just a small 
amount of material can thus be lost during this phase (as 
M oc M at this phase). This clearly prevents any large 
-•^^F yield. With the new lower M values, the same star 
enters the WC phase with a mass larger than 43 Mq , and 
more than 27 Mq are lost during the WC phase. This 
favours higher ^^F yields. 

6. Conclusions 

Revised ^^F yields from non-rotating WR stars and the 
first evaluation of the yields from rotating such stars are 
presented. The new yields in absence of rotation are sig- 
nificantly lower than those of Paper I, as illustrated by 
the Z = 0.02 60 Mq case, where the reduction amounts 
to more than a factor of 30. Rotation does not help in this 
matter, and even reduces the yields. This drastic decrease 
of the predicted ^^F yields mainly results from the adop- 
tion of reduced mass loss rates, and to a lesser extent from 
the selection of a smaller core overshoot. 

Taken at face, these new predictions discard WR stars 
as important sources of the galactic ^^F. Let us however 
emphasise that they suffer from uncertainties originating 
from at least two sources. As discussed above, p™'^ is very 
sensitive to (1) the still poorly known ^^F (a , p) ^^Ne rate 
(see also Stancliffe et al. l2005|l . (2) the still uncertain wind 
mass loss rates. Interestingly, Meynet & Maeder (2005J 
suggest on the basis of arguments concerning the WR pop- 
ulation at solar and higher than solar metallicities, that 
the mass loss rates during the post-core H-burning WNL 
phase might be underestimated. Higher mass loss rates 
during this short stage would uncover more rapidly the 
He-core and would likely favour the ejection of ^^F. 

All in all, we consider that the question of the con- 
tribution of WR stars to the galactic "'^^F remains largely 
open. It appears reasonable at this point to refrain from 



drawing any far-reaching conclusion based on the present 
WR yield predictions, particularly in attempts to build 
galactic chemical evolution models. 
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